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Abstract 

In this paper we study the fusion frame potential, that is a generahzation of the Benedetto- 
Fickus (vectorial) frame potential to the finite-dimensional fusion frame setting. The structure 
of local and global minimizers of this potential is studied, when we restrict the frame potential 
to suitable sets of fusion frames. These minimizers are related to tight fusion frames as in the 
classical vector frame case. Still, tight fusion frames arc not as frequent as tight frames; indeed 
we show that there are choices of parameters involved in fusion frames for which no tight fusion 
frame can exist. Thus, we exhibit necessary and sufficient conditions for the existence of tight 
fusion frames with prescribed parameters, involving the so-called Horn-Klyachko's compatibility 
inequalities. The second part of the work is devoted to the study of the minimization of the 
fusion frame potential on a fixed sequence of subspaces, varying the sequence of weights. We 
related this problem to the index of the Hadamard product by positive matrices and use it to 
give different characterizations of these minima. 
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A Hadamard products and indexes. 
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1 Introduction 

Fusion frames were introduced by P. G. Casazza and G. Kutyniok under the name of "frame of 
subspaces" in [8]. They are a generahzation of the usual frames of vectors for a Hilbert space Ti; 
indeed frames of vectors can be treated as "one dimensional fusion frames". During the last years, 
the theory of fusion frames has been a fast-growing area. Several applications of fusion frames have 
been studied, for example, sensor networks [H], neurology [l9], coding theory [1], [5] , [E], among 
others. Specifically, applications which require distributed processing can be well described and 
studied using fusion frames. We refer the reader to the work by P.G. Casazza, G. Kutyniok and 
S. Li [To] and the references therein, for a detailed treatment of the fusion frame theory. Further 
developments can be found in [7], [9] and [20] . 

In finite dimension, a fusion frame is a sequence of subspaces of F" (F = C or M) together with a 
set of positive weights such that the weighted sum of the orthogonal projections to these subspaces 
(called fusion frame operator) is a positive invertible operator (see Definition 12. l.ip . As in the case 
of vector frames, it is usually desired this invertible operator to be a multiple of the identity. In 
this case, the fame is called a tight fusion frame, also noted TFF. However, tight fusion frames 
might not exist for a fixed choice of dimensions d = {di)i^i^ of the subspaces, for any sequence of 
weights w = {wi)i^i^^^ (see the discussion following Proposition 13.1.1] ). 

Inspired on the work in classical frame theory (see [21 dS])) we study a convex functional on 
fusion frame operators (the FF potential or FFP), also studied by Casazza and Fickus in [7], 
which generalizes the Benedetto-Fickus frame potential. In this paper we analyze its local and 
global minima. This is motivated by the fact that local (global) minimizers characterize unit norm 
tight vector frames (see [21 U] [S] [U US]). Since the FF potential can be seen as a "measure of 
orthogonality" of the frame vectors, it provides also an interesting geometrical description of fusion 
frames. These considerations motivate the study this type of minimizations in the fusion frame 
context. A related study can be found in [7], a work that the authors became aware in an advanced 
stage of writing this paper. 

The main tool used in [TH] for these problems for vector frames, namely majorization of matrices, 
can be replaced in the context of fusion frames by the theory developed by Horn and Klyachko 
in order to have a spectral characterization of hermitian matrices which are the sum of a set of 
hermitian matrices. For example, this approach provides necessary and sufficient conditions of 
existence of TFF's, summarized in a family of inequalities. This technique, although seems to be 
rather impractical due to the complex conditions involved, becomes an useful tool in the study of 
the spectral structure of the FF potential minimizers. 

We first consider the problem of existence of TFF's. We show some dimensional restrictions 
on the subspaces regarding this problem, and then give equivalent conditions for the case of fixed 
dimensions and weights. We refer to [7] for further developments in this direction. 

The rest of the paper deals with the minimization of the FF potential on some sets of Fusion 
Bessel sequences (i.e. sets of projections and weights whose fusion frame operator is not necessary 
invertible). Mainly, we work on Bessel sequences with fusion frame operator of trace one. This is 
a natural restriction in order to avoid scalar multiplications, and it allows an interpretation of the 
FF potential as a measure of how far is the fusion frame operator to the suitable multiple of the 
identity corresponding to the (possibly non-existing) tight fusion frames with trace one. A detailed 
discussion of this approach can be found in subsection 12.21 
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The minimization of the FFP is done in three different settings: first, by fixing the weights w 
and the dimensions d of the subspaces. Then, fixing only the dimensions d. Finahy, we consider a 
fixed sequence of subspaces {Wi)i^i^, and optimize over the set of admissible weights. In the three 
cases, the minimization is made under the previously mentioned "trace one" restriction. 

For the first problem, a geometrical approach similar to that done in [18] allows us to obtain a 
characterization of local minimizers of the FF potential: they are orthogonal sums of tight frames 
on each eigenspace of the frame operator. Then, using Horn and Klyachko techniques, we prove 
that all minimizers (even those which are local minimizers) have the same eigenvalues, with the 
same multiplicities. Similar results are obtained for the second problem (fixing only the dimensions 
of the subspaces). 

The last section of the paper is devoted to the study of the optimization of the fusion frame 
potential of fusion frames obtained from a fixed sequence of subspaces in C"" which generate the 
whole space. Since every sequence of weights makes it a fusion frame, we seek for the best choice of 
weights, meaning those which minimize (globally) the FF potential. Then, we establish a connection 
between this optimization problem and the Hadamard indexes studied in [T2] (which involve the 
Hadamard or entry-wise product of matrices) of a kind of Gramm matrix associated to the fixed 
subspaces. Using these tools, we get a characterization of the set of optimal weights, and a way 
to compute them under some reasonably assumptions on the initial sequence of subspaces. This 
analysis seems to be new even for the case of vector frames. 

However, as it is shown by an example, the minimizers could "erase" some of the initial subspaces 
(i.e. the set of optimal weights could have zeros). Moreover, it is possible to obtain a minimizer 
which is a Bessel sequence of subspaces which is not generating, a phenomenon which does not 
happen in the previous settings. Motivated by this problem, we study the geometry of the set of 
all weights w which minimize the FF potential, and in particular their possible supports (namely, 
those sub-indexes i such that > 0). At the end of the section, we present some examples which 
illustrate these type of anomalies. 

The paper is organized as follows: Section 2 contains preliminary definitions on fusion frames 
and the basic notation used throughout the paper. This section also contains a brief exposition 
of Horn-Klyachko's compatibility inequalities. Section 3 is devoted to the study of minimizers of 
the FF potential, restricted to the sets of fusion frames detailed before. In Section 4 we analyze 
the problem of minimize the FF potential for a fixed sequence of subspaces, varying the weights. 
The paper ends with an appendix containing definitions and several results concerning Hadamard 
indexes of positive matrices, which are related to the contents of Section 4. 

We wish to acknowledge to P. Casazza and M. Fickus for letting us know about their excellent 
work [7], which is closely related with the present paper. 

2 Preliminaries and Notations. 

In this paper A^„(C) denotes the algebra of complex nxn matrices, Ql{n) the group of all invertible 
elements of A^„(C), U{n) the group of unitary matrices, A^„(C)sa (resp. Mn{^)ah) denotes the 
real subspace of hermitian (resp. anti-hermitian) matrices, A4„(C)"*" the set of positive semidefinite 
matrices, and gi{n)+ = Mn{'C)+ ngi{n). Given T G A^n(C), R{T) denotes the image of T, N{T) 
the null space of T, o"(T) the spectrum of T, trT the trace of T, and rk T the rank of T. 

Given m G N we denote by Im = {1, . . . ,rfi} and 1 = Im G M™ denotes the vector with all its 
entries equal to 1. If f € M", we denote by diag(7;) G A^„(C) the diagonal matrix with v in its 
diagonal, and by G M" the vector obtained by re-arrangement of the coordinates of v in non- 
increasing order. If T G A4„(C)sa , we denote by A(r) G M"- the vector of eigenvalues of T, counted 



3 



with multiplicities, in such a way that A(T) = X(T)K 

Given a subspace W C C", we denote by P\y g ^An{C)~^ the orthogonal projection onto W, i.e. 
R{Pw) = W and N{Pw) = W'^. For vectors on C" we shall use the euclidean norm, but for 
matrices T G A4n{C), we shall use both the spectral norm ||T|| = ||T||sp = max ||Tx||, and the 

||x-|j=l 

Frobenius norm ||T||2 = (tr T*T)^ — { ^- This norm is induced by the inner product 

{A, B) = tr B*A, for A,B Mn{C). 

2.1 Frames of subspaces, or fusion frames for C" 

We begin by defining the basic notions of fusion frame theory in the finite dimensional context. 
For an introduction to fusion frames for general Hilbert spaces, see [8], [10] or |20j . Briefly, a 
fusion frame for is a generating sequence of subspaces, equipped with weights assigned to each 
subspace. Nevertheless, we prefer to give the "frame style" definition, which adjusts better to our 
purposes. 

Definition 2.1.1. Let W = {VFjjgi^ be closed subspaces of 7i = C", and w = {tUjjiGi™ G K>o- 
The sequence Ww = {wi , W^i)iei„i is a fusion frame (FF) for H, if there exist A,B>0 such that 

A\\ff < ^ wf \\Pwjf < B\\ff for every / G W . (1) 

If only the right-hand side inequality in ([T]) holds, then we say that Ww is a Bessel sequence of 
subspaces (BSS) for 7i. The frame operator of W«, is defined by the formula 

Sw^ = Y.wlPw,eMn{C)+ . (2) 

Observe that Ww is a FF if and only if S'yy^ G Ql{n)'^ and, in this case AIn < 5'w„ < B In ■ We 
say that is a tight FF (TFF) A = B , in other words, if Syv^ = AIn ■ A 

Remark 2.1.2. Let W^, = (wi , Wi)iei„ be a BSS for C". For each i G I m , we can take an o.n.b. 
Bi = {e^j^}j^j^ of Wi . Hence, for every / G 7Y, we have that 

PwJ=Y.(f'4^)^? ' foriGl^ =^ Sy,J=Y,^lPwJ=Y.Y.^f^'^^^f)'^^^f ■ 

jeJi i&m. i&mj&Ji 

Therefore, Wu, induces a vector Bessel sequence T = {wi e^-^ : i ^ Im , j ^ Ji} which has a very 
useful property: Its frame operator 5jf = Sy\>^ . A 

2.2 Sets of fusion frames and the FF-potential 

We shall establish several notations regarding sets of FF's and BSS's : 

Notations 2.2.1. Fix n, m G N and consider a Hilbert space H = C^. 

1. We shall denote by Sm,n the set of all FF's of the form Ww = {w, W) = {wi, Wj)jg][^, where 
w G M>o ^'^d ^ ^ generating sequence of subspaces of TC. 
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2. Given a sequence d G N™ such that trd = ^ dj > n, we denote by 

Sm,nid) = {Ww & Sm,n ■ dim Wj = di for every i G 1^} , (3) 
Similarly, we denote by Bm,n{d.) the set of BSS's with the same dimensional restrictions. 

3. Given v G M™o , we denote by 

Bm,nid,v) = {Ww ^ Sm,n{^) ■■ W = v} and Sm,n{d,v) = Bm,n{d,v) nSm,n{d) , (4) 

the subsets of Sm,n(d) and Sm,n{d.) with a fixed sequence of weights v. 

4. Finally, we denote by 

Bl^M) = G e™,n(d) : trSw» = E = l} • (5) 

and 54,„(d) = i3^,„(d) n 5„,„(d) . 

5. We say that a pair (d, w) G N'" x ]R™q is normalized if tr d > n and E wf di = 1. Observe 
that (d,-u;) is normalized if and only if Bm,n{^-,w) ^ 'Sm.n(^)- ^ 

The following definition is suggested by the classical Benedetto-Fickus potential, whose value in a 
vector frame J- can be calculated as FP(jr) = tr S"^ . 

Definition 2.2.2. Given a ESS W«, = {wi^Wi)ii=j^ , the Benedetto-Fickus fusion frame potential 
(FF-potential) of is given by: 

m 

FFP (W^) = w^w] tiiPw.Pw,) = tr S^^ . (6) 

Notice that in this case FFP (Ww) = FP(J^), for any vector Bessel sequence J- obtained from Ww 
as in Remark 12.1. 2 [ We define also the following matrix: 

m m 

Pg(W») = Yl "'^'^i " ^ ^ >fn(C)+ . (7) 

i,j=l i,j=l 

The matrix Pg(-) is related to the so-called q-potential [17] defined in the more general context of 
reconstruction systems. Notice that the Benedetto-Fickus fusion frame potential can be computed 
in terms of Pg(W^), since FFP (W^) = trPg(W^) . A 

The scope of this paper is to study minimizers of the FF-potential. In order to avoid scalar 
multiplications (note that FFP (Wt-w) = ^^FFP(>V^^,)) we shall restrict ourselves to minimize 
the FF-potential on subsets of 5^ „(d) or B^^ „,(d), for d and m fixed. In other words, we shall 
minimize the FF-potential for those frames Ww such that tr S'w^ = 1. 

This specific restriction is justified because, if there exist tight FF's in 5m.n(d), then their FF- 
potential and frame bounds are determined exactly by the trace of their frame operators. Namely, 
if Ww € 5^,n(d) is tight, and tr 5w„ = a, then 5w„ = ^ In and FFP (Ww) = ^ ■ 

Even in the case that there are no TFF's in 5m,n(d), this restriction seems to be quite natural. 
Indeed, for BSS's with fixed trace, the FF-potential can be seen as a measure of the (Frobenius) 
distance of their frame operators to a fixed multiple of the identity: 
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Proposition 2.2.3. LetWw e5^,^„(d). Then 

\\- In-S^^^\\l = tTS^^ - - = FFPiW^) - - . 
" n 112 » ^ ^ 

Proof. Since tr S'w^ = 1 , a direct computation shows that 

\\— In - Sw^Wl = tr ( ^ In •S'w^, + 5^ ) = tr - — . □ 



The last result shows that if there exist tight FF's in 5^ n(^)) then they are the unique global 
minimizers of the FF-potential on 5^ n{'^)- ™ the case that there are no TFF's in 5^ nl^)' the 
minimization of the FF-potential becomes more interesting: it provides the elements of 5^ n(d) 
that can be expected to have the best properties. 

In this paper we deal mostly with these type of minimizations under two different further restric- 
tions: we work in the set Bm,n{^iw) Q n(d) for a fixed normalized pair (d , w), or we fix a 
generating sequence W of subspaces, and minimize the FF-potential over all sequences w G M>q 
such that Ww € nCd). 

2.3 Klyachko-Fulton approach 

Recall that given x € M", we denote by € M" the vector obtained by re-arrangement of the 
coordinates of x in non-increasing order. Given x, y M" we say that x is subniajorized by y, and 

k k n n 

write X y, if ^ < ^ yj^ for every fc S I„ . If we further have that tr(x) ■= Yl — Vi 

i=l i=l i=l i=l 

then we say that x is majorized by y, and write x ~< y. 

Example 2.3.1. As an elementary example, that we shall use in what follows, let x € M>q and 

< a < tr(x) < b. The reader can easily verify that ^ 1„ x bei . A 

(Sub)majorization between vectors is extended by T. Ando in [1] to (sub)majorization between 
self-adjoint matrices as follows : given A, B G A^„(C)sa , we say that A is submajorized by B, and 
write A -Ky^ B, if A(^) ~<w \{B). If we further have that tr(A) = tr(i?) then we say that A is 
majorized by B and write A ~< B. 

Although simple, submajorization plays a central role in optimization problems with respect to 
convex functionals and unitarily invariant norms, as the following result shows (for a detailed 
account in majorization see Bhatia's book [3]). 

Theorem 2.3.2. Let A, B A^„(C)'"^. Then, the following statements are equivalent: 

1. A B. 

2. For every unitarily invariant norm || • || in A1„(C) we have ||A|| < \\B\\. 

3. For every increasing convex function / : M ^ M we have tr f{A) < ti f{B). 

Moreover, if A -<„, B and there exists an increasing strictly convex function / : M ^ R such that 
tr f{A) = tr f{B) then there exists U G U{n) such that A = U*BU. □ 
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In what follows we describe the basic facts about the spectral characterization of the sums of 
hermitian matrices obtained by Klyachko [K] and Pulton [T3]. Let 

C = {Ul, • • • , jr) G : jl < i2 • • • < jr} ■ 

For J = (jl, . . . ,jr) E /C", define the associated partition 

X{J) = {jr-r,...,h-l). 

Denote by LR^{m) the set of (m + l)-tuples (Jq, . . . , Jm) G (/C")'""''^, such that the Littlewood- 
Richardson coefficient of the associated partitions A(Jo), . . . , X{Jm) is positive, i.e. one can generate 
the Young diagram of A( Jq) from those of A(Ji ),..., A(Jm) according to the Littlewood- Richardson 
rule (see [13]). With these notations and terminologies we have 

Theorem 2.3.3 ([EKll]). Let A^ = Aj = (A^, . . . , A^*^) G M" for i = 0, . . . , m. Then, the following 
statements are equivalent: 

1. There exists Ai € Mn{C)sa with X{Ai) = Aj for < i < m and such that 

Ao = Ai + ... + Am . 

2. For each r G {1, . . . , n} and ( Jq, . . . , Jm) ^ LR^{m) we have 

m 

E4°'sEE^f (8) 

plus the condition J2 "^j J2 Yl "^j ■ 

j=l i=lj=l 

Moreover, if (Aj)^o ^^'^ item 1. above and (Jq, . . . , Jm) £ LRp{m) satisfy equality in 

then there exists a subspace L C C" with dim L = r, that simultaneously reduces Ai for < i < m 
and such that X{Pl Ai) = (Aj*'*)jg j^, where Pl denotes the orthogonal projection of C" onto L. □ 

We shall refer to the inequalities in ([8]) as Horn-Klyachko's compatibility inequalities. 

3 On the existence of tight fusion frames. 

The following facts exemplify the difference between the theory of vector frames and that of frames 
of subspaces. In [2] (see also [6] and [18]) it is shown that the local minimizers of the frame potential 
on the set 

Fi,n = {^= {fih&m ■■ each /, G C" and tr(S'^) = H/.f = 1 } 

are tight frames. Since the set „ is compact and the frame potential is a continuous function, 
there must be global (and hence local) minima of the frame potential. This was used to give an 
indirect proof of the existence of such frames in the vectorial case. 
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3.1 Dimensional restrictions 



Let d G N"^ with trd > n and consider the set n(d) defined in Eq. ([5]). Using Remark 12. f .21 it 
fohows that if d = 1^ , then we can identify ni^) with „, and the previous comments can 
be apphed. Hence it seems natural to ask whether there always exist TFF's in n(d)> since they 
would be all the global minimizers of the FF-potential on ni^)- The following results show that 
in general the answer is no. 

Proposition 3.1.1. Let (d, E N"* x be a normalized pair, with M = trd > n, and assume 
that G 5^ „(d) is a TFF, so that Syv^ = ^ In ■ If there exists i G Im such that 

M — di = dk < n — 1 =^ wf = — and Pwj = for every j G Im \ {«} • 

Proof. Consider the tight vector frame T = {wf ej*^ : i G Im , J £ } associated to Ww , as 
described in Remark I2.1.2[ Let G G Mi\j{C)'^ denote the Gramm matrix of the vector frame J- 
and let Gi = wf I^. denote the Gramm matrix of each subsequence {wf e^*^ : i G 1^.}. Then each 
Gi is a X di principal sub-matrix of G. By Cauchy's interlacing principle [3] we get: 

Aj(G) > \j{Gi) > XM^d,+j{G) for 1 < j < d^ , 

where A(G) = (Aj(G))jgn^^ (resp A(Gj) G W^^) denotes the vector of eigenvalues of G (resp. Gi) 
counting multiplicities and with its entries arranged in non-increasing order. By assumption, 

Aj(G) = ^ for j Gin , Aj(G) = for n<j<M, and Xj{Gi)=wf for j Gld, ■ 

Thus ifJ2dk<n-l, then M - + 1 < n and ^ = AA./-d,+i(G) < Xi{Gi) = wf < Xi{G) = ^ . It 

(i) 

is known that, in this case, each of the vectors e^ , 1 < j < di , must be orthogonal to every other 
vector in the -tight vector frame which implies the last assertion of the theorem. □ 

Example 3.1.2 (About the existence of tight frames in nCd)). Consider now d = (2,2) and 
assume that there exists Ww G iS2 3(d) that is tight. That is, we assume that there exist two 
subspaces Wi C with dim Wi = 2, i = 1,2 and wi , W2 £ M>o such that ^ I^ = wf Pwi Pw2 ■ 
Since di , ^2 < 3 — 1 we conclude from Proposition 13.1.11 that wf = W2 = , and Pwi Pw2 = 0; 
which is impossible. This argument can be extended to show that if the choices of d = {di)i^i^ 
are such that each di is relatively small compared with n and ^ dk then there are no TFF's in 

5^ „(d). For example, in 5^^ 2^-1(2 • Ifc), 5| 7(3,3,3), 5| 9(4,4,4), etc, there are no TFF's. A 

Remark 3.1.3. The previous results show some dimensional restrictions for the existence of TFF's 
in n(d). In the paper by Casazza and Fickus [7], some sufficient conditions on n, m and d are 
given (particularly if d is a multiple of Im); which assure the existence in 5^ n(d) of such fusion 
frames. For further results in this direction, see also [16]. A 

3.2 Characterizations for fixed weights 

The following theorem gives us some general bounds for Pg(-) and states several conditions on 
Ww G 'Bm,n(d) which are equivalent to the assertion that Ww is a -^--TFE. 



8 



Theorem 3.2.1. Let Ww G Bm,„(d), where '^i ^ 1- ^^^n 

^/^^Pg(W»). (9) 
For every u.i.n. \\ ■ \\ on A4n{C) with associated symmetric gauge function ip we have that 

^ < l|P,(>v»)ll ■ (10) 

For every increasing convex function f : ]R>o M with f{0) = we have 

n-/(^) <tr /(P,(W^)) . (11) 

Finally, the following conditions are equivalent: 

1. is a^-TFF. 

2. Majorization holds in 

3. There exists u.i.n. \\ ■ \\ such that equality holds in (jlOp 



4- There exists an increasing strictly convex function f : ]R>o — > M>o with /(O) = such that 
equality holds in (fTTj) . 



Proof. Since tr (S'wu,) = J2 wfdi>l then (see Example 12.3. ip it follows that \ln <w Sy\!^ and 
hence 

- = tr (i Inf < tr S^^ = tr Pg(W^) =^ \ln Pg(W^ ) . (12) 

Notice that by Theorem 12.3.21 then (llOp and (llip are consequences of this last fact. Assume that 
majorization holds in Q, so then we have 

tr In = tr Pg(W^) = tr 5^^ . 

Since 5w„ and the function f{x) = x"^ is strictly convex, by Theorem 12.3.21 we conclude 

that there exists a unitary U € U{n) such that = U*{^In)U = ^In- On the other hand, if 
there exists an u.i.n. || • || such that equality holds in (fTOll then, using the right-hand side of (fT2]) 
we get 

^ ^(1) = ||P,(W.)|I > ^LZiiZM ^(1) > J_ ^,(1) , (13) 

which implies that tr (-^ /„)^ = tr Pq(yV«,). As before, we conclude that 5'w„ = -i- /„ . Similarly, if 
there exists an increasing strictly convex function / : ]R>o M>o with /(O) = such that equality 
holds in ()11|) then Theorem 12.3.21 and the right-hand side of (I12p imply that Sw^, = — • Finally, 
it is clear that in case Ww is a TFF then Pg(Wu,) = In ■ The last part of the theorem follows 
from this fact. □ 

Theorem 3.2.2. Let (d, w) € N™ x ]R™q be a normalized pair. Then, the following statements are 
equivalent: 

1. There exists Ww = {wi , Wi)i^i^ S 5^ „(d) which is a ^ - TFF. 
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2. For every 1 < r < n — 1 and every (Jq, . . . , Jm) € LR^[m) we have that 

- < V u;M Jin{l,...,(ii}|. 
n ^-^ 

3. There exists an orthogonal projection P € Alm(A4n(C) ) with tr(P) = n and such that, if 

= {wi Wj Pij)i,jei„, with Pij £ Mn{C) for i,j £ , then 

—5— Pa = ( —5— Pa] = ( Pa ] and tr ( Pa ] = di for every i Gim ■ 
wf n n J \wjn J \wf n J 

Proof. Notice that the condition in ([T]) is equivalent to the existence of orthogonal projections 
{Pi}i£i^ such that tr(Pj) = dj for 1 < i < m and such that ^ wf Pi = ^ /. Hence, by Theo- 

rem 12.3.31 it follows that condition ^ should hold, since these are Horn-Klyachko's compatibility 
inequalities for the spectra of {vjf Pi}i^i^ and ^ . The converse of the previous implication also 
follows from Theorem 12.3.31 since a self adjoint A operator with A(^) = (a, . . . , a, 0, . . . , 0) € is 
necessarily of the form A = aP for some projection P £ A^„(C). 

Assume now ([T]) and let Pi = Pwi for 1 < i < m, where Ww = {wi,Wi)i^i^ is a ^-TFF. 
Let us consider Vi G A4„(C) a partial isometry such that V*Vi = Pi, 1 < i < m. Define V* = 
[wiV{\w2V^\ ■ ■ ■ \wmV^] G Mn,m-n{C) and notice that V*V = ^ ^fy^Vi = i^- Hence VV* = 

{■WiWjViV*)^j^^ G MmiMniC)) is such that VV* = ^ P ioi an orthogonal projection P = 
iPij)Tj ^ ■Mmi-Mni'C) ). By comparing the diagonal blocks we get that 

- Pa = wf ViV* =^ Pu = ViV* , for every i £lm ■ 

n WJ n 

Conversely, if P = (Pij) G A^m(7V4„(C) ) is as in ^ then there exist matrices Vi G A^n(C) such that 
if V* = [wiVil ■ ■ ■ IwmVm] G A^n,m-n(C) then P = VV* (since rank P = n). But then, comparing 
the block diagonal entries we get that V^Vi is an orthogonal projection with tr(-^ V*Vi) = di for 
1 < i < m. Hence V*V = In , that is 

Y.wf^V*V = -In . 

Therefore, if we let Wi be the range of V* we get that Pwi = ^V'*Vi so then W = {wi, Wj)jgn^ is 
a J- -TFF for C". □ 

n 

4 Minimization for fixed weights 
4.1 Lower bound for the potential 

In this subsection we translate, using Remark 12.1.21 some well known results about vector frames 
(see [6] or [20]) to the FF context. An interesting fact is that there is a notion of irregularity, 
defined in terms of the parameters of a given FF, which agree with the vectorial n-irregularity of 
their associated vector frames. Nevertheless, the lower bound obtained for the FF-potential is not 
always attained in the set I3m,ni^,'w) (see Example I4.1.3P . 



10 



Definition 4.1.1. Given a pair d € N™ and w = £ '^>o ' consider its q-irregularity defined as 

j m 

Jo(d, w) = max ^ j e : (n - ^ dj) w| > ^ 
if this set is not empty, otherwise Jo(d, w) = 0. 

Proposition 4.1.2. Lei (d, G x R™q 6e a normalized pair, where w = wK Recall that 

m 

M = di>n. Let jo := Joid,w) and c = < w'j^ . IfWw G Bm,n{d,w) then 

i=l 

jo m 

FFP{W^) >Y^d^wf + {n- ^ di) c^. (14) 

i=l j=io + i 

Moreover, equality holds in (|14p i/ and onZy i/ t/ie following two conditions hold: 

1. Pw, Pwj = for 1 < i ^ j < jo and 

2. {wi , Wj^^.^+i is a TFFfor span{Wi : 1 < i < jq}^. 

Proof. 1. Let Ww £ Sm,n(d, and let ^ = {wf €>^^ : i € Im > j ^ 1^^} be an associated vector 
frame, as described in Remark I2.1.2I Let a € M^^^ denote the vector whose coordinates are the 
norms of the elements of JF arranged in non-increasing order. Then, 

k-l 

Wk = o,M(k,j) , where M{k,j) ='^di+j , for k elm and j € Id^ ■ 

1=1 

We now consider the n- irregularity r„(a) of the vector a : 

M 

r„(a) =max|j G : {n - j)aj > ^ } , 

i=j+i 

jo 

if the set on the right is non empty, and r„(a) = otherwise. It is straightforward that r„(a) = di 

i=l 

in the first case. Therefore, inequality ()14p can be deduced from [6l Theorem 10] (see also |20j). 
The same result of [H] shows that equality in Eq. (fM|) implies that 5i = {e^j^ : 1 < « < io ) J ^ Id, } 

is an orthonormal system in C", and S2 = {e^*^ '■ jo + ^ ^ i ^ , j £ 1^, } is a tight frame for S-^, 
where S = span {Si} = span {Wi : 1 < i < jo}. □ 



The following example shows that the lower bound in ()14p is not sharp in general. 

Example 4.1.3. If we set n = 3, d = (2,2) and wi = W2 = ^ then ^ wf di = 1 and Jo{d,w) = 0. 

Therefore, by Theorem 14.1.21 the equality (I14p holds only in tight fusion frames. Still, there are no 
tight FS in S'2,3(d,fi;) C S'2,3(d) since the Example 13.1.21 shows that there are no tight FS in the 
(bigger set) 52,3 (d). A 
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4.2 Structure of local minima: The geometrical approach 

In what follows we consider a perturbation result for Bessel sequences of subspaces. We begin 
by considering some well known facts from differential geometry that we shall need below. In 
what follows we consider the unitary group U{n) together with its natural differential geometric 
(Lie) structure. It is well known that the tangent space 7}^ U{n) at the identity can be naturally 
identified with the real vector space 

Mn{C)ah = i ■ Mn{C)sa = {X G Mn{C) : X* = -X} , 

of anti-hermitian matrices. Given G G A1„(C)"^ we consider the smooth map 

^^G-l^in) ^U{G) CMn{C) given by "ifGiU) = U*GU , U ell{n) , (15) 

where hl{G) is the unitary orbit of G. Under the previous identification, the differential of at 
a the point /„ £ U{n) in the direction given by X G A^„(C)ah is given by 

{D^g)i„ {X) =XG-GX = [X, G] . (16) 

It is well known that the map ^g is a submersion of U{n) onto U{G). Therefore, the differential 
{D'^g)i„ is an epimorphism, and hence (fT6l) gives us a description of the tangent space of the 
manifold U{G) at the point G: We have that TgU{G) = { [X, G] : X G Mn{C)ah }• 

Let us fix some notations. We denote by 

Mi{C)sa = {A& Mn{C)sa : tr A = 1} and Ml{C)sa = {Ae Mn{C)sa : tr A = 0} . 

Observe that 7W^(C) sa is an affine manifold contained in the real vector space M-niC) sa ■, whose 
tangent space is the subspace M.^{C)sa ■ On the other hand, given X,Y ^ M.n{C)sa , it is easy to 
see that tr XY G M. Therefore the inner product {A, B) = ti B* A of 7W„(C) still works as a real 
inner product on Ain{'C)sa ■ 

Given a set {Pj : j G Im} ^ ■Mn{C)sa of projections, we denote by 

{Pj ■■ j G ImY = {Ae MniC) : APj = A for every j G I^} . (17) 

Note that {Pj : j G Im}' is a closed selfadjoint subalgebra of 7Wn(C). Therefore, the algebra 
{Pj : j G Im}' ^Clfi <J=^ there exists a non-trivial orthogonal projection Q G {Pj : j G Im}'- 

Theorem 4.2.1. Let Ww = {wi , Wi)iei,^ G B}^ ni^^'^) B'^^- Denote Pj = Pwj for every 
J G Im • Let ^ : U{n)^ — > M\{C)sa ^ A^n(C)sa be the smooth function given by 

m m 

^{Ui,...,Um) = Yl ^] UjPj Uj = w] {Uj) , for (f/i , . . . , f/m) e U{nr • 
j=i i=i 

Then the following conditions are equivalent: 

1. The differential of ^ at I := (In, • • • , In) ^ U{n)^ is surjective 

2. {Pj :i Gim}' = C/n. 

m 

In this case, the image of contains an open neighborhood of'^{I) = ^ '^'i Pj ^ni'^)sa, o,rid 
^ admits smooth (and hence continuous) local cross sections around ^{I)- 
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Proof. It is clear from its definition that ^' is a smooth function. Moreover, under the previous 
identification TiJJ{n) = M.n{C)ah > and using Eq. (fTOj) . we can see that 

m 

D^j{X,,. ..,Xm)=Y^ w] {Xj,P.^ , for (Xi, ...,X„,)e Mn{C)2 ■ (18) 
j=i 

The tangent space of A4l^{C)sa is the real vector space A4^{C)sa , which has a natural inner product 
given by {Y, Z) = tr(yZ). Denote by T = D^j and assume that T is not surjective. Then there 
exists 7^ y G ^A^{C)sa which is orthogonal to the image of T. Using Eq. we deduce that, 
for every {Xi, . . . , Xm) G A4„(C)J^ , it holds that 

m m 

i=i i=i 

Since each [Pj , F] G A4n(C)ah , we can choose each Xj = [Pj ,y], and so Eq. (fT9]) implies that 
[Pj , y] = for every j G Im • In other words, that Y G {Pj : j G Im}'- On the other hand, since 
7^ y G MniC) sa and tr y = 0, then Y ^ CIn ■ The converse follows from the previous argument, 
by taking Y G {Pj : j G Im}' such that 7^ y = y* and tr y = 0. □ 

Let {d,w) G x M'"q be a normalized pair. We shall consider on Bm,n{<i,'w) the distance 

dpiWy, ,WJ= max \\Pwi - Pw' II 

(recall that the weights are fixed), called punctual, and the pseudo-distance 

called operatorial. The problem of finding local minimizers for the FF-potential can be stated for 
anyone of those distances between EE's. 

Corollary 4.2.2. Let {d,w) G x W^q be a normalized pair. Assume that Ww G ^^.^(d,^;) 
satisfies that {Pwj '-3 £ Im}' = Cln ■ Then is a FF and the map 

S : ]3m,n{d,w) Ml{C)sa given by S'(V^) = •S'v^ = ^ WiPVi , 

i&m. 

for Vw = (ujj, Vi)jgi,„ G Bm.,ni(i,w), satisfies that 

1. The image of S contains an open neighborhood of Sy\!^ in A4}^(C)sa ■ 

2. S has dp -continuous local cross sections around Syy^ . 

Proof. Notice that the condition {Pwj : J G Im}' = C/„ implies that W is a generating sequence 
of subspaces. To prove the properties of the map S", just compose a local cross section for the map 
^ of Theorem 14.2.11 (which is open in /) with the map <1> : hl{n)^ — > n.(d) given by 

$(C/i , . . . , Um) = {wi , Ui{Wi) )i6n„ = [w, , R{UiPw,U*) . 

Observe that S o = , so that S is open in $(/) = Ww ■ □ 
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Remark 4.2.3. Here is an alternative statement of Corollary I4.2.2t Under the same assumptions 
and notations about , it holds that 5w„ S Gl{n)~^ and, for every sequence {Sk)keN ■M.li{C)sa 
such that Sk > Sy\>^ , there exists a sequence {Vk)k>ko iii '5m, n(d , w) such that 

dp{Vk , ) > and Sx;,^ = Sk for every ko < k . 

k—*oo 

This formulation of Corollary 14.2.21 generalizes \18\ Thm 5.3] to the context of fusion frames with 
fixed weights. A 

It is not clear that a dp -local minimizer for the FF-potential on Bm,n{d, w) must be a fusion frame, 
i.e. its frame operator is an invertible operator. The following Lemma shows that this is true. 

Lemma 4.2.4. Let {d,w) e N™ X MJ^Q be a normalized pair. Let Wto be a dp-local minimizer for 
the FF-potential in Bm,n{*i,'w). Then S'wu, invertible {equivalently, Wto € Sm^ni^^w) ). 

Proof. Suppose that 5w^ has nontrivial nullspace N{Sy^^^ ). If x S N{Sy^^ ) 

= {Sw^ x,x) = ^ {w]Pwj ^^x) =^ \\Pwj x\\^ . (20) 

j&m j&m 

In other words, Wi C R[Sw^ ) for every i € Im • Since tr(d) > n > dimi?(5w„ ), we deduce that 
there exists i 7^ j in 1^ such that Pwj Pwi 7^ 0. Fix that pair i, j. Fix also / G VF?, \ Wj- and 
g € N{Sy^^ ) two unit vectors. For every t € [0, vr/2], take the unit vector g(t) = cost • / + sint • g. 

Let Vi = WiQ span {/}, Wi{t) = © span {g{t)} and (t) € Bm,n{d, w) the sequence obtained 
by replacing Wi by Wi(t) in Ww . As 5 € Wj^ for every A; € Im , for every t G (0, 7r/2], 



^ '— ^ ^ = }2 wfwi I tr(Pw'^ PwJ - iT{Pw,{t)PWk) 

k=l ^ 
m / 

= E ^Mi^^Pw. Pw,) - tr {g{t)g{tr Pw,) - tr{Pv,Pw,) 
k=i ^ 

m / 

= E wWk ^<Pw. Pw,) - cos2 t tr (//* Pw,) - tr{PvM 
k=i ^ 

> EwM[^<Pw.Pw,)-tr{ff*Pw,)-tr{Pv.Pw,))=0 



k=l 

because tv{ff*Pwj) = \\Pw, ff + 0. Hence FFP(V^(t)) < FFP (W^ ) for every t € (0,7r/2]. 
Taking t ^ 0, we have that VV«, (t) -5 , and this contradicts the minimality of W^u . □ 

Given S € M.n{S^)sa with (j(S') = {^1 , . . . , /x,.}, we denote by P^^y^iS) = PN(S-iikin) ^ A^n(C)+, 
the spectral projection of S relative to fi^ 5 for A; G 1^ . These projections satisfy that 

p 

1. Pfj_^ (S) Pfj,. (S) = if fc / j, and ^ P^^. (5) = (i.e., they are a system of projectors). 

k=l 

P 

2. For every fc G 1^ , it holds that 5P^^(5) = fik P^j,(5), so that S = Yl f^k P^lki'^)■ 

k=l 

The following theorem generalizes a similar result given in the paper by Casazza and Fickus [71 
Theorem 4], for the case w = Im ■ Nevertheless, our approach is based on completely different 
techniques. 
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Theorem 4.2.5. Let {d,w) S N™ x W^q be a normalized pair. Let Ww G Bm,n{d.,w) be a local 
minimizer of the FF-potential with respect to the distance dp . If a{Sy^^ ) = {/ii , . . . , Hr}, then 

P^^k(Sw^ ) e Cw„ = {Pwj ■■ 3 G Im}' for every keir . (21) 
The same property holds whenever is a dp-local minimizer in ni^)- 

Proof. Recall from Lemma [4.2.4l tliat € 5m, n(d, fu), in other words that ^ (t(S'w^ ). Consider 
the set Q of finite systems of projectors {Qk}k&p such that each Qk € . Observe that Q is 
not empty because {In} € Q. Then Q has a maximal element {Qfcjfceip with respect to the order 
induced by refinement. Fix /c G Ip . For each i G Im put Mi = Wj n R{Qk) and A/i = Wi n R{Qk)'^ . 
Using that Qk G Cyy^ , we get that each Wj = Aii © A/i . Set = dim A4j and r = (ri , . . . , rm). 
Then, the sequence Wk,w = (w^i > ■^j)iGim is a FF for R{Qk). We claim that Wk,w is a local 
minimizer of the FF-potential in 

B{Qk,r,w) := [Vu, = (wi , Vi)i<^i^ G ^^^^^(r,^;) : Vi C R{Qk) for every i G I,„ } . 

Indeed, given Vw G B{Qk,'r,w), put = (■Wi , ffiA/"j)j6i^ G ,Bm,n(d,u;). Observe that the map 
Vw ^ Vw preserves the distance dp . Moreover, since Qk G Cy\>^ , then each Pj\f. = {In — Qk)PWi 
so that, by Eq Q, FFP (V^) = FFP (V^) + FFP ( {wi,Ni)i(zi^ ), and the second summand does 
not depend on Vw ■ Then, the claim follows from the fact that Wk,w = • 

Observe that S'w,,, commutes with Qk ■ We now show that Sy\!^ Qk = Ok Qk for some Uk G cr{Sy\>^ ) . 
Indeed, by the maximality of {Qi}^^i in Q, it follows that there is no non-trivial sub-projection 
Q' of Qk such that Q' G {Pmj '■ j G Im}'- Then we can apply Corollary 14.2.21 (taking TL = R{Qk) 
and renormalizing the traces) to show that every positive operator (with the correct trace) near 
Swvj Qk has the form 5v,„ for some G B{Qk,r,w) close to Wfc . But if 5w„ Qk is not a scalar 
multiple of Qk , then we can choose S\;^ in such a way that 

FFP (V^) = tr S^^ < tr ( S^^Qk ) = FFP (W,.) . 

But this contradicts the fact that is a local minimizer of the FF-potential in B{Qk, r, lu). Hence 

p 

Qk = ctk Qk and Qk < Pa^iSw^ ) ■ Using that ^ = I„ , it is easy to see that each 

k=l 

Pi^,{Sw^ )=^Qk^ Cyv^ , where Ji = {k elp : ak = □ 

fce Ji 

Remark 4.2.6. Next we give two reinterpretations of Theorem 14.2.51 Under the assumptions and 
notations of the theorem, the following properties hold: 

1. For each i G Im , there exists an o.n.b. Bi = {e^*'' : j G 1^-} of Wi , consisting of eigenvectors 
of 5w„ • Indeed, observe that each P\y^ = Ylk&i PWiPfiki'^Ww ) and the fact that, for a fixed 
i G Im ) the projections Pwi Pfiti^Ww ) are pairwise orthogonal. 

2. For each /i^ G cr{Sy^^ ), denote by M.k,i = ^{Sww ~ l^k In) H Wi . Then, it holds that the 
sequence Wk = {wi , ■Mk,i)i&Tn is a tight FF for N[Sy\;^ — fik In)- This follows because its 
frame operator SyVk = -P^fc i^w,, ) Sw^u = l^k P^ik ('S'w„ )• A 
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4.3 The eigenvalues of all ds-minimizers coincide 

Recall that, given S G Al„(C)sa , we denote by X{S) € M" the vector of the n eigenvalues of S, 
counted with multiplicities, in such a way that A(«S') = X{S)^. 

Lemma 4.3.1. Let {d,w) € N™ x M™q be a normalized pair. Then 

1. Bm,n{d.,w) is dp-compact. 

2. The set Am,n(d, w) = {A(5w„) : G Bm^ni^, w)} is a convex and compact subset of M". 

3. The set {Syy^ : € Bm,nid,w)} is compact and closed under unitary conjugation. 

Proof. Let W^, G Bm,nid,w)- For t € [0, 1] consider A = t A(5>v„) + A(5n)^) and notice 

that A = A-l-. Therefore, for every admissible (m + l)-tuple (Jq, . . . , Jm) £ LR^{m), 1 < r < n — 1 
we have 

■m m 

since both A(S'wu,) and \{Sy\;iJ) satisfy Horn-Klyachko's inequalities. Hence, by Theorem 12.3.31 
there exists Vw = {wi, yi)i&m € •Sm,n(d, w) such that A(S'v„) = A. This shows that A € Am_„(d, if), 
so that Am,n((l, It;) is convex. The fact that the set {S'w : Ww S 'Sm,n(d, w)} is closed under unitary 
conjugation is apparent. Finally Bm,n{d^w) is dp-compact because each Grassmann manifold 

VdAn) = {P = P*P e Mn{C)+ ■.tTP = d^} = {UP,U* : U € U{n)}=U{P{) , 

for every fixed Pi E Vdi{n) , is compact. This follows because U{n) is compact. By continuity, the 
other sets involved are also compact. □ 

Theorem 4.3.2. Let {d^w) € x M™q be a normalized pair. Then, 

1. The spectra (with multiplicities) of all the frame operators of global minimizers of the FF- 
potential in Bm,n{d,w) coincide. 

2. The local minimizers of the FF-potential in Bm,n{d-, w) with respect to the pseudo-distance ds 
lie in 5m,^„(d, if) and are also global minimizers. 

Proof. Let Ww G ^rn.,n(d, w;) and notice that FFF (Ww) = ||A(S'yv;^) |p. Since Am,n{d.,w) is a 
compact convex set, then there exists a unique Aq S Am^„(d, u;) which minimizes the euclidean 
norm on A^, n(d, w). Hence, if Ww is a global minimizer of the FF-potential in Bm^ni<i, w), we can 
conclude that ||A(5w„)|P < ||Ao|P, which implies that A(5w„) = Aq . 

Observe that the map a : Bm,n{d,w) — > Am,^„(d,t(;) given by cr{Ww) = A(5w„) is continuous 
with respect to the pseudo-distance ds of Bm^n{d,w). Moreover, a is an open map. Indeed, 
fix Ww G Bm,nid.,w), A = X{Sy\>^) G A^, n(d, and take G A^, n(d,tf) close to A. Let 
U G U{n) such that Syy^ = i/diag (A ) C/*. By Lemma f4.3.H there exists Vw G Bm,n{d,w) such 
that Sv^ = J7diag(/i) U* . Now observe that ds{Ww ,Vw) = \\Swu, - Sv^W = ||A - ^||oo • 

Therefore, if is a d^-local minimizer of the FF-potential in Bm,n{d,w), then A(S'w^,) is a local 
minimum for the euclidean norm in the set A^, n(d, w) Q M". By a standard computation, the 
convexity of A^, n(d,u;) implies that A(S'w^) must be the global minimizer Aq , and therefore Ww 
is a global minimizer in Bm,n{d, w). □ 
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Conjecture 4.3.3. Local minimizers of the frame potential in Sm,.n(d,i«) (resp. n(d)) with 
respect to punctual distance dp are also global minimizers. 

In some particular cases (i.e. for particular choices of the parameters n, m, d and w), there is an 
affirmative answer for this conjecture (see [TJ Theorem 5]). 

Remark 4.3.4. All the previous results remain true if one replaces Bm,n{d.,w) by „(d) (i.e., 
minimizing without fixing the sequence of weights). We present some of the new statements without 
proofs, since they are based on techniques that are similar to those already developed. 

1. As in the proof Lemma Horn-Klyachko's compatibility inequalities ([8]) show that 

A;'„ „(d) = {X{Sw^) '■ G ^m,ni^)} is compact and convex. 

2. This fact implies that the (ordered) spectra of the frame operators of global minimizers of 
the FF-potential in Bj^ „(d) coincide. 

3. Finally, the argument of the proof of Theorem 14.3.21 can be adapted to yield that local 
minimizers of the FF-potential in B^^ n(d), with respect to the operator distance ds , are also 
global. A 

As in the case of fixed weights (Lemma l4.2.4p . the global minimizers for the FF-potential in n(d) 
are fusion frames. 

Proposition 4.3.5. Let Ww o, ds-local (and hence global) minimizer for the FF-potential in 
n(d)- Then its frame operator 5w„ G £//(n)+. In other words, Ww & n(d)- 

Proof. Let J = {i Im '■ Wi ^ 0}, and k = ^ di . Note that, if A; > n, we can apply Lemma 14.2.41 

(fixing the weight wj) and we are done. 

We assume that k < n, and will obtain a contradiction. Without loss of generality, we can suppose 
that J = Ir . It follows immediately that Wi _L Wj for 1 < i ^ j < r and wf = ^ for every 
l<i<r. Thus, FFP (Wu, ) = 1 ■ Moreover, if d = dr+i , then k + d > n. Otherwise, if we take a 

subspace Wr+i ^ ( |J Wi) with dim Wr+i = d and we set w'^ — k+d ' ^ ^ -Ir+i ! then we get 
a BSS in „(d) with FF-potential -j^^ < ^ . 

Therefore, we can construct V„ = {vj , V^)jei,.+i G B^^ „(d) in the following way: 



Wj if 1 < i < r 

^~\b ifi = r + l ~ U ©W^i]^ e r ifj = r+l 




V? = ■( and Vi 



1=1 



where T C ©[^^^ Wi is a subspace with dimT = d + k — n (so that diml4-_|_i = d), and ka + db = 1. 
It is easy to see, by taking an orthonormal basis of each subspace Vi , that 

f{a) = FFP (V„) = ka^ + db'^ + 2{d + k- n)ab with ka + db=l , a G [ , ^ ] . (22) 

2(n — k) 

Easy computations show that /'(-|) = ^ > 0. Since f{\) = \ , there exist pairs (a, b) such 

that the FF-potential is lower that , which contradicts the minimality of yV«, . □ 
Next, we summarize the facts described in Remark 14.3.41 and Proposition 14.3.51 

Theorem 4.3.6. Let d E N"'. Then, local minimizers of the frame potential in n(d) with 
respect to the pseudo-distance ds lie in „(d) and are also global minimizers. Moreover, the 
spectra of the frame operators of these local minimizers coincide. □ 
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5 Minimization for fixed subspaces 



In this section we shall characterize the sequences of weights which minimize the potential of a fixed 
sequence of subspaces. The main tools are some results about Hadamard indexes of [12], which 
we shall state in some detail in the Appendix. Recall that, for A,B£ A^„(C), their Hadamard 
product is the matrix Ao B = (A^ Bij)ij^i^ £ ^An{C). 

5.1 Minimal weights 

In this section we fix jtt, € N, a Hilbert space Ti. with dim 7^ = n, and d G with M = tr d > n. 
We also fix a sequence W = {VFijjgi^ of subspaces which spans H, such that each dim Wj = di . 
Our aim is to minimize the FF-potential over all sequences w G such that Ww G n(^)- 

Recall that the Benedetto-Fickus FF-potential of Ww is given by: 

m 

FFP (W^) = ^ wfw]ti{Pw,Pw,) = tiS^^ . (23) 

Definition 5.1.1. Let Ww = {wi , W^Oieim ^ ^rn,n{^)- We denote by 

B = Bv\!^ G A^m(]R) the matrix given by Bij = wjvJj ti^PwiPwj) ■ A 

Lemma 5.1.2. Let >V„, = {wi , Wi)iei„ be a FF for H. Then, By^^ G 7W„(R)+. 

Proof. By^^ G A^m(IK)^ because it is the Gram matrix for the vectors {wfPwi}i&,n ™ the euclidean 
space Aln(C) with the inner product defined as {X , y) = tr Y*X. □ 

Notations 5.1.3. We shall fix some notations and assumptions: 

1. We begin with a fixed normalized sequence of weights, in the sense that 

-1 

w = {wi)i(zi^ is given by Wi = d-'^ for every i G Im • 

-1 

Observe that the condition Wi = d • ^ means that each "vector" Wi Pwi of Wyj has size 
llwj -PiyJIa = 1- This justifies the word "normalized" of . 

2. Given a sequence of weights a = (aj)ien^ G M™, we denote by a • the Bessel sequence of 
subspaces a • Ww = {ai Wi Wj)igi„ . 

3. In a G M!p , then tr(5'a.w„ ) = X] ^I'^l'^i ■ Therefore, as we start with normalized weights, 

a • G 5i,„(d) ^ Yl = = 1 ^ ll^ll = 1 • (24) 

i 1_ 

4. Let A = Ayv^ G A1m(M) the matrix given by Aij = )ij = WiWj tr{PwiPWj) , for 
hj £ Im ■ Observe that A is selfadjoint, but possibly A ^ J^m(M)~^- On the other hand, 

m 

FFP (Wu,) = E iT{PwTw,) = Pw„ 11^ and also 

FFP (a • ) = E a2 w] t^{Pw.Pw, ) = ||aa* o Ayy^ \\l . 
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5. Using the previous identities, we can now define the main notion of this section: 

/(>V„, ) := min FFP (a • ) = min ||aa* o Ays>^ \\l . A 

\\a\\=l \\a\\=l 

In order to compute I{Ww ) as well as to describe the set of weights a S I^>o; with ||a|| = 1 for 
which ) = FFP (a • Ww ), the main tools are some results about Hadamard indexes of [12], 

which we shall state in some detail in the Appendix. Here we just give the basic definitions. 

Definition 5.1.4. Let B G Mm{'C)~^ . The minimal-Hadamard index of G is 

7(5) =max { A > : 5oC > AC for every C G X™(C)+} . 

For A G Mra{C) sa ) we define the spectral and the || • Hj Hadamard indexes: 

Isp{A) = min ||Aoxx*||sp and I2{A) = min ||^°2;x*||2 . 
11x11=1 l!3^ll=i 

For a matrix G G A1„(C), we write ^ G if all entries Gij > 0. Given J C I„ with \ J\ = k we 

ij 

denote by Gj G Alfc(C) the submatrix of G given by Gj = (Gij)..^j. Similarly, if x G C", we 
write X ^ if X G M!p and xj = {xj)j,zj G C^'. A 

From the previous definitions and Eq. ([251) . we get the fundamental equality: /(W«, ) = /2(^w„)- 
Now it is clear why the results of the Appendix can be useful for computing I(yV«,). 

Remark 5.1.5. Let a G M™. Then 

m 

FFP (a • ) = ^ a^i a] wfw] tT{PiPj) = {{B o aa*) a, a) = (S(a o a*), (a o a*)) (26) 

where, as before, B = By^^ G M.mi^)'^ ■ Moreover, by Eq. ([32]) in Proposition IA.1.2t 

I{y\^w) = min oaa*||^ = min ||Boaa*|| = Isp{B) . (27) 

||a||=l ||a||=l 

This identity is useful because ^ G A1„i(M)^ and its spectral index is easier to compute. Indeed, 

ij 

observe that tr P]v^ Py/j > for every i,j G Im , since each Pwi G ^An{C)~^. A 
5.2 Critical Points and local minimizers. 

In what follows, we shall use all the assumptions and notations of the previous subsection, but we 
need the following extra notations: 

1. Given a G W", we write z = a o a = {af , ■ ■ ■ , a^) and J = supp{a} = {j G : 7^ 0}. 

2. Let L : ^ M+ be given by L(a) = FFP (a • ) = \\A^^ o aa*||^ = z , z). 

3. We consider the affine manifold Aq = {x G M'" : trx = 0} and the compact convex simplex 
A = {x G Ao : x^O}. 

4. S"'-^ = {x G : ||x|| = 1} is the unit sphere of M™. 

Lemma 5.2.1. Let O^a G S'^~^. Denote B = Byv^ ■ The following conditions are equivalent: 
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1. a is a critical point for the map L restricted to the sphere 5"™ ^ . 

2. B jzj = I{B j)l J , where J = suppja} . 

m 

Proof. Observe that i^(a) = {Bz, z) = ^ij '^1 • Since B = B* = B^, we have that 

m m 

= VL(a) = 4 (ai ^ hij a| , . . . , a„ ^ 6^^ a^) = 4 o a . 

The tangent space of S""~^ at a is {a}-*". Hence, a is a critical point for S""~^ if and only if 

= (VL(a), y) = 4 {Bz o a, y) for every y G {a}"*" <^=> Bz o a G span {a} . 
This is clearly equivalent to the equation Bjzj = Alj , for some A G M. In this case, since O^B 

ij 

with < Bii for every i G and O^z, we can conclude that A > 0. Moreover, by Proposition 
IA.1.21 applied to the matrix Bj , we have that 

Bjzj = Xlj =^ ^J^ = '^J ^ /(i?j)-i=tr ^ = (lj, ^) = A-Ml, ^) = A^^ 

Therefore A = I{Bj) and Bjzj = I{Bj)lj . □ 

Theorem 5.2.2. Let O^a G S"^~^ such that Bjzj = I{Bj)lj , i.e., a. is a critical point of L 
restricted to 5™"^. Then, the following conditions are equivalent: 

1. I{Bj) = Isp{B) = I{W^). 

2. a is a global minimum of L restricted to 5™"^. 

3. a is a local minimum of L restricted to S"^~^ . 

4- Bz ^ I{Bj)l. In other words, that {Bz)j > I{Bj) for every j ^ J. 

Proof Denote A = Ay^^ . Recall that hiAf = Isp{B), by Eq. ([27]). By Lemma EXIl 

\\Aoaa*\\l = {Bz,z) = I{Bj){l,z)= I{Bj) . 

This gives the equivalence 1^2. Observe that, if b G /S™"^, then t/; = bob G A. For each G A, 
consider the line 7^ : [0, 1] — > A joining z and w, given by the formula 7u;(t) = (1 — t)z + tw, for 
every t G [0, 1]. Consider the map pw : [0, 1] M"^ given by 

p^it) = {B^n,it) , 7n.it)) = (1 - t)^ {Bz ,z)+t^ {Bw , w) + 2t(l - t) {Bz , w) , (28) 

for every t G [0, 1]. Since {Bz , z) = I{Bj), the derivative of p^ evaluated at zero is 

p,^(0) = -2 {Bz , z) + 2 {Bz , w) = 2(^{Bz , w) - I{Bj) 
On the other hand, for every t G M, 

p^,{t) = 2{Bz, z) + 2 {Bw ,w)-4:{Bz,w) = 2 {B{z - w) , z - w) > . (29) 
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Since is a second degree polynomial, its leading coefficient is ^ Pw{t) ^ 0. Suppose now that 
Bz ^ I{Bj)1l. Using that it; G A, we get that 

{Bz , w) > I{Bj) (1 , w) = I{Bj) =^ p^(0) > ^ p^{t) > for every t > . 

Therefore Pw{^) > Pw{0)- In other words, we have proved that {Bw , w) > {Bz , z) for every w & A. 
This implies that a is a global minimum of L restricted to S""~^. 

Suppose now that a G S"^~^, z = a o a, and that there exists A; € Im such that {Bz)k < I{Bj). 
Observe that a 7^ , because {Bek)k = bkk = -^(-^{A;})- 

Let w = e/fc G A, and consider the curves 7^ and pw defined before. By the previous computations, 
we have that Piu(O) = 2( {Bz)k — I{Bj) ) < 0. Therefore, for every t > small enough, we have 
that 7w(i) S A and (Bj^it) , 7w{t)) < {Bz, z). Taking the vectors a(t) = sgn(a) 7,u,(t)^ G 5™"-*^, 
we conclude that a fails to be a local minimum. □ 

Remark 5.2.3. Given a G S"^'^ , then z = a o a G A and L(a) = {Bz, z). Therefore 

Uy\)^) = 12(A) = min ||Aoaa*||a= min ||yl o aa* lln = min (Sz, z) , 

since every z G A produces a unit vector ^ a = G Then in order to get the unit vectors 

a which attain this minimum, it suffices to characterize the sets 

SiyVw ) = argminKi^z , z)} and J(W^ ) = ^ Im '■ J = supp{z} for some z G 5(W^ )} . 

If Im ^ J{y^w ), it is possible to obtain minimizers a • >V«, which are not fusion frames, because 
5a.w„ i Gl{n)+ (see Example EZZD- Still, if /(W^ ) < then 5a.w„ G ^/(n)+ for any 

minimizer a, since in such case Proposition 12.2.5] implies that ||/ — nSg^.yy^ \\ < 1. Otherwise, the 
characterization of the set J{yVw ) is useful in order to discern if there are optimal sequences of 
weights a such that a • Ww remains being a FF. Item 4 of Theorem 15.2.21 gives a description of the 
elements of 5(W^ ). But its proof gives more information: A 

Corollary 5.2.4. Consider Ww & Sm.n{d.), A = A-^^^ and B = Byy^ as before. Then 

1. The set 5(>V.u, ) = argmin{(i3z , z)} is convex. Moreover, for any point zq G 5(W^ ), 

S{W^ ) = {zo + N{B) ) n A . 

2. J{yVw ) is closed under taking unions, so that Jy\!^, = IJ J{Ww ) = U suppjz} is an 

element of J{Wui), CLnd there exists z\ G 5(yV«, ) with maximal support. 

Proof. 1. Let z, if G A, and consider the function, defined in Eq. ()28l) : 

Pz,wit) = {B{ (1 - t)z + tvj) , {l-t)z + tvj) , t G M . 

Suppose now that w,z G 5(W«, ) and w z. Using that pz^^ is of second degree, the equality 
'Pz,w{f) = 2 {B{z — vu) , z — w) > given by Eq. ([2^]) . and the fact Pz,w{t) > for every 
i G M, we can conclude that 

Pz^w is constant 'pz,w{t) = <^=^ z — w G N{B) . 

On the other hand, we have that Pz, wi^) = Pz,wiO) = niin^ Pz,w{t). This implies that the 

map pz,w is constant, so that ^{t) G S{y\^w ) for every t G [0, 1], and z — w G ^(B). The 
proof of the fact that (zq + N{B) ) n A C S{Wnj ) for every zq G S{Ww ) is similar. 
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2. Let z and w in iS(W«,), with supports J\ and J2 respectively. Then, since the entire hne 
+ (1 — G 5(>Vtt, ) (t G [0, 1]), if we take u = + (1 — l)w for any t € (0, 1), it is easy to 
see that u G 5(W«, ) and suppjti} = Ji U J2 . Since J(yV^ ) is finite, also the set 

•^w„ = U supp{z} G J(W^ ) . 

ze-s(w„) 

Hence Jyv^, is the support of some z\ G 5(yV«, ). □ 
Corollary 5.2.5. Let B = Byy^ and A = Ayy^ as before. Assume that there exists v G such 
that u ^ and B v = \. Denote a = (tr v) 2 [vi , . . . , ). Then 

||a|| = 1 and (tr w)"^ = FFP (a • >V„, ) = ^2(A)^ = ) . (30) 

Proof. The fact that u ^ and Bv = \ implies, by Propositions IA.L2] and lA.1.31 that 

h{Af = hp{B) = I{B) = I{Bj) = {iiv)-^ , 

where J = supp{i)} = supp{a}. Since z = G A, then a. = ^ gm~i ^^^^ _ 

Hence, by Theorem [5221 we have that z G 5(W^ ) and a G S'"'"^ satisfies Eq. ([30]). □ 

Remark 5.2.6. The results of this section seems to be unknown still for the case of vector frames. 

In this case our restrictions translate to the following: Let J- = {fi)ieim ^ frame for Ti. such 

-1 

that each ||/j|| = 1 (i.e., di = 1 =^ Wi = d^^ = 1). For a G M"^, we consider the sequence 
Sl ■ T = {ai fi)i£i^ , and we define = min FP (a • JF). Then, all the results of the section 

l|a|| = l 

remain true if one consider the matrices 

A:F=(\{fj,f,)\) eMmiC)sa and i?^ = ( |(/,. , ) g >I™(C)+ . 

Some proofs are slightly easier in this case, because 1{T) = l2{_Ajr) = l2{Gj^), where Gjf is the 
Gramm matrix of T: Gjp = ( {fj , fi) )^ G A4m{C-)~^. Observe that the diagonal entries of the 
three matrices involved are equal to 1. A 

/4 1 3\ 

Example 5.2.7. Let B = i 1 4 2 . Since A = {B^^.)ij(z^ G ^/(3)+, we deduce that A is 

\ 3 2 4 / 

the Gramm matrix of a Riesz basis J- of C'^. Let u = ( | , | , 0) ^ 0. Observe that 

Bv =13 =^ (by Corollary [523]) z = {trvy^ v = (^l ,1 ,0^ e S{J=') . 

Since N{B) = {0}, Corollary 15 . 2 .41 assures that iS(^) = {z}, and Jq = {1, 2} is the maximal support 
for S{J-). Taking a = z^ , we have that a • ^ is the unique scaled sequence of J- with minimal 
Frame Potential, but it fails to be a frame for C^, because it has just two non zero elements. A 

Example 5.2.8. It can be proved that every G G A^3(C)^ such that rk G = 2 and Ga = 1 for 
every i G I3 (considered as the Gramm matrix of a frame J- for with three unitary elements), 
satisfies that the minimizers a • JT of the BF-potential are frames for C^. 

Indeed, given z G S{J-), it is easy to see that J = suppjz} has more than one element (otherwise 
z = ei for some z G I3). If J = I3 there is nothing to prove. Assume that suppj^} = J with 

I J| =2. If rk Gj < 2 we must have Bj = f | | j . In this case, I{Bj) = Isp{B) = 1. But the 
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unique matrix B S M^IC)'^ which satisfies that O^B, Isp{B) = 1 and Ba = 1 for every i S I3 is 



ij 



B = \- t*. Indeed, if some By < 1, then Isp{B) < hp ^ ^ = < 1. Finally, since 

1 = rk 1 • 1* = rk i3 = rk G o G > rk G = 2, we have a contradiction. A 

Remark 5.2.9. Let W = {W^ijieim be a generating set of subspaces. Given a partition {Jk}k&p 
of the set 1^ , we say that the sequence {y^k}k&p of W given by Wk = iWi)ieJk ^ partition in 
orthogonal components of W (briefly POC) if Wi _L Wj for every pair i £ Jk , j £ ^m \ Jk ■ 

Note that by definition the trivial partition given by Ji = Im produces a POC of W. If {yVk}keip 
is a POC of W, we say that it is maximal if the only POC of each Wfc is the trivial one. It is clear 
that there always exits such a maximal POC for W. 

Let {Wfcjfegip be a maximal POC of W with | J^j = rrifc for 1 < A; < p. Let G ]R™''= be such that 
||afc|| = 1 and /(Wfc) = FFP (a^ • Wfc) for each 1 < k < p. Then, there exists 7 = {'yk)keip S 
with II7II = 1 and such that 

p 
k=i 

Conversely, if a = (ai, . . . , ap) with ^ a^ G M™fe and ||a|| = 1 is such that I{W) = FFP (a • W) 
then afc ^ and /(Wfc) = FFP (||afc||~^ a^ • Wk), for 1 < A; < p. Hence, we can restrict our study of 
the optimal weight of sequence of subspaces to each of the components of the maximal partition. 
This in turn implies that we can reduce the problem of describing the optimal weights to the case 
where the matrix B (which has non-negative entries and is positive semi-definite) is irreducible i.e., 
none of its symmetric permutations can be written as the direct sum of two matrices. This last 
property is relevant in the theory of matrices with non- negative entries. 



A Hadamard products and indexes. 

In this section we recall some definitions and results from [T^] which are closely related with the 
problems of Section 4. The exposition is done with some detail for several reasons: a) Most results 
we state are explicitly used in the previous section, b) The formulation of these results given in 
|12j is quite technical and intricate, so we intend here to give a clarified version, c) Although 
some results in the appendix are not directly applied, they are included since they give effective 
criteria for computing the indexes and the vectors that realize them. This is relevant since we have 
identified these objects as the optimal weights and the minimal potential for fusion frames. 

A.l Basic definitions and properties 

We begin with an extended version of Definition 15.1.41 
Definition A.1.1. Let G € Mm{C)+. 

1. The minimal-Hadamard index of G is the number 

/(G) =max{A > : GoB > AS for every B £ Mm{C)+} ■ 

2. Given an u.i.n N in A4m{C), the A^-Hadamard index of G is 

InIG) = max {A > : N{G o B) > \N{B) for every B G 7W„,(C)+} 
= min {N{G o B) : B £ 7W„(C)+ and N{B) = l] . 
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The index of G associated with the spectral norm || • || — || • is denoted by Igp{G^^ and the one 
associated with the Frobenius norm || • ||2 is denoted by l2{G). A 

Proposition A. 1.2. Let G G Alm(C)+, 1 = (1, 1, . . . , 1) G and E = 1 • 1^. 

1. I{G) 7^ if and only if 1 G R{G). If there exists y G C'' such that Gy = 1, then 
I{Gr^ = E = (y, 1) = PiG^m = mill { {Gz, z) : ^ z, = I } . 

/ \-l det(7 

If G > 0, then also /(G) = ( ^ {G^^' ^ 



det(G + E) -detG 

2. /(G) < In{G) for every u.i.n. N. 

3. If J C I„ , then /(Gj) > /(G) and In{Gj) > In{G). 

4. If D = diag (d) G A4m(C)+ is diagonal, then iNiD) = N'{d~'^y^. In particular, 

m = h,{D) = ( ^ dr'Y' and h{D) = ( ^ df) ^ . 

5. Both indexes I2 e Isp are attained by matrices B G A4,n(C)''' of rank one. This means that 

hiG) = min ||Goxx*||2 and Isp{G) = min ||Go2/2/*|| . 

11^11=1 Il3/ll=l 

Moreover, the minima are also attained at vectors x ^ (or y^O). 

6. Let 5 = GoG G >(,n(M)+. Then /a (G) = /.p(S)i = /^^(GoG)^ . 

7. Moreover, if ^ S G 7Wm(M)+ and ^ G Mn{^)sa is given by = 5^/^ for I < i, j < m 

ij 

then, even if A ^ A^m(I^)^5 the index /2(^) of Definition 15.1.41 still satisfies 

12(A) = min ||^oxx*|L = min ||/?oxx*||^ . (32) 
||x|i=i" |lxi|=i" 

8. It holds that Isp{G) = inf { Isp{D) : G < D and D is diagonal } . Therefore 

/2(G) =ini { {J2d^^) ' : < /) is diagonal and G o G < } . □ 

1 

Proposition A. 1.3. Let G G A4m(M)+ such that G. Then /sp(G) = /(G) ^ there 

ij 

exists n ^ such that Gn = 1 . □ 

Proposition A. 1.4. Let G G A4,„(C)+. Denote by P = G o G. UxeR^, then 

II G o xx* II 2 = IGjjPlxjplxj 1^ = (P(x o x), X o x) = ((P o XX* )x, x) < HP o xx*|| . 
id 

Take x^O such that ||x|| = 1 and ||Goxx*||2 = /2(G). Then 

(G o G o XX* )x = I{Pj)x , where J = {i G Im : Xi ^ 0} . 
In this case, they hold that 

1. The vector u = I{Pj)~^{x j o xj) G C"' has strict positive entries and Pju = Ij . 

2. I,p{P) = hp{Pj) = I{Pj). 

3. Isp{P) = ||Poxx*||. □ 
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